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We combine linear response theory and dimensional regularization in order to derive the dy- 
namical Casimir force in the low frequency regime. We consider two parallel plates moving along 
the normal direction in D— dimensional space. We assume the free-space values for the mass of 
each plate to be known, and obtain finite, separation-dependent mass corrections resulting from the 
combined effect of the two plates. The global mass correction is proportional to the static Casimir 
energy, in agreement with Einstein's law of equivalence between mass and energy for stressed rigid 
bodies. 
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^ I The Casimir force on moving plates has nontrivial properties, which have been studied in recent years. It provides 

■ for dissipation of the plate's mechanical energy, an effect closely related to the emission of photon pairs g] g. This 
coupling also gives rise to decoherence of the plate's motion As far as the photon emission effect is concerned, it 
is interesting to assume that the plate oscillates with a frequency close to a resonance frequency of a cavity resonator, 
because the effect is greatly enhanced in this case Q]. However, to achieve such unusually high mechanical frequencies 
is a challenge not yet undertaken. More generally, dissipative effects involve field modes with frequencies of the order 
of the mechanical frequency, and are very small when the latter is well below the lowest cavity frequency. 

On the other hand, dispersive effects result from the contribution of the entire field spectrum, and hence are less 

■ sensitive to the variation of the mechanical frequency. In this paper, we consider the Casimir forces on two parallel 
moving plates, in the nonrelativistic and long-wavelength approximations. Photon production rates for this geometry 

. have been derived in this limiting case ||] . Here we focus on the dispersive force, and take an additional approximation: 
' we assume the motion to be slow in the time scale associated to the time-of- flight of light between the plates. In this 
■~»>^ , limit, the main nontrivial corrections to the static Casimir force are inertial forces, i.e., contributions proportional 
' to the accelerations of the two plates. This effect was analyzed in detail by Jaekel and Reynaud ^ in the case 
' , of a one-dimensional space. Here we consider D spatial dimensions, and use dimensional regularization, as in the 

formalism employed by Ambjorn and Wolfram for the static Casimir effect 0. 
(-H ' The Casimir force and photon emission rates for a moving dielectric semi-infinite half-space (single interface) in D 
• • . dimensions were also computed ^ . Generalizing the result of a previous work , Ref. ||] obtained cut-off dependent 
mass corrections for the moving interface, which renormalize the bare mass. The corrections and the bare mass are 
both infinite (in the no cut-off limit), and when considered separately, unaccessible to measurements, only the total 
\^ mass possessing a physical meaning. The mass corrections we shall derive here have a different interpretation. Rather 
5^ , than mass renormalization, they represent a measurable physical effect, not present in the case of a single plane 
interface (or plate) considered in |^ The corrections are finite, and depend on the distance between the plates. 
As required for consistency, they vanish at the limit of large distance, leaving the place for the unperturbed free- 
space masses of each plate. The free-space mass eventually also incorporates the vacuum radiation pressure coupling 
through a renormalization of a bare mass, in the line of Refs. |^] and Here we skip this important problem to focus 
on the effect brought about by the presence of the second plate, and start from the experimentally known free-space 
masses The mass of a single plate may also be modified when its surface is corrugated. Mass corrections for one 
and two parallel corrugated plates were computed in the case of lateral motion . 

This paper is organized as follows. In section 2, we derive most of the formal results (leaving technical aspects to 
three appendices), and some physical implications are discussed in section 3. Additional comments and conclusion 
are presented in section 4. 
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II. INERTIAL FORCES 



We consider a massless scalar field cj) in Z?— dimensional space. The field obeys Dirichlet boundary conditions at 
the moving plates xi — Sqi(t) and xi = a + Sq2{t) : 

cj){&ii{t),X2,...,XD,t) = = (j){a + Sq2{t),X2,...,XD,t) 

Hence Sqi (t) and Sq2 (t) represent the time-dependent deviations of the positions of the plates from their unperturbed 
values, which correspond to a Z3-dimensional static cavity of length a. 

The vacuum radiation pressure force on plate a (a = 1,2 are the labels for the two boundaries) is written as 

fa^ + Sfa(t), where fa'^ is the static Casimir force operator for the unperturbed cavity length a. We use hats to 
distinguish between the force operators and their average values. For instance, the static average forces satisfy 
/i"^ = if 2°^) = — but no such simple relation holds for the force operators themselves. We compute Sfa(t) up to 
first order in Sqi{t) and Sq2{t) and their time derivatives, from the fluctuations of the force in the unperturbed case, 
with the help of linear response theory. In the frequency domain (we employ capital letters to denote the Fourier 
transforms), we have 

2 

SFa{uj) = ^Xat3{uj)Qf3{uj), (1) 

13=1 

with the Fourier transform defined as: 



The susceptibility Xafii^) may be computed with the help of linear response theory. In this approach, the response is 
derived from the fluctuations of the force in the static (unperturbed) case . The imaginary part of the susceptibility, 
which accounts for dissipation, is given by 

ImXa/3(cj) = ^ (Ca/3(t^) - Ca/3(-w)) , (2) 

where Ca[}{Lo) is the Fourier transform of the correlation function of the force: 

The average is taken over the vacuum state, and for the unperturbed configuration corresponding to the stationary 
cavity of length a. The force is computed from the stress tensor (we take c = 1 in this section) 



Stj = -di4>dj4> + 



.k=l 



(3) 



with i,j = 1, ...,£) 
When considerii 
stationary configuration, and hence vanishes at = and Xi 



When considering the fluctuations of fa ^ we replace (p by the unperturbed field (po, which corresponds to the 



(poixi,. ..,XD,t) ^y^Yl \ ^ tD-1 [ — L j exp(ifc{„^} • a;||)a,.{„^} exp(-zcj;^{„^.}i) +II.C., (4) 

1=1 {„a V \ a J 

where x\\ = {x2^ ...,a;£)), and {uk} = {n2, is a list of integer numbers (we have taken periodic conditions at 

the boundaries of a cell of measure L^^^ on the hyper-planes at = and xi = a; as usual the limit L — *■ oo 
is assumed), and the corresponding sum is over all integer values. Each list {uk} corresponds to a wavevector 
^{"fc} ~ 27r(n2, ...,njj)/L parallel to the hyper-planes. Hence a given field mode corresponds to an integer I and a list 

{rife}, and its associated frequency is cui^^nk} — \J Q-^/o-Y + ^{nj.}- 
The force on the boundary at xi = is given by 

fT\t)= I d^-ix||5i4(0+,:r||) = -i / (a,,(/.o(0+,x,|,t))% (5) 

J xi=0 ^ Jxi=0 
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whereas the force on the second boundary is 



/fW-^/ d^-S (a,,^o(a-,X||,t))', (6) 

^ Jxi=a 

the integrals taken over the quantization ceU. The spectra of force fluctuations are computed in Appendix A, yielding, 
with the help of (H), the imaginary parts of the four susceptibilities. They satisfy the symmetry relations xii(w) = 
X22(w), Xi2(t^) = X2i{^)- We find 

ll — l ^2 — 1 

x[S{uj- (uji, + Lui,)) -5iuj + {lui, + uji,))], (8) 
where sii(/i,^2) = 1, s^ihjh) = and 

Lui = ^{In/aY + kl (9) 

According to Eq. (|^), for a given mechanical frequency w, dissipation comes from contributions of cavity modes 
satisfying the resonance condition associated with the emission of photon pairs: mi^ +toi2 — ^- In fact, the mechanical 
power dissipated by the force equals the total radiated power, as expected by energy conservation 

On the other hand, the real part of the susceptibility, which accounts for the dispersive component of the force, is 
computed from the dispersion relation (P J denotes Cauchy's principal value) 

ti,eXaf3[UJ) = -V 1 duj ; , (10) 



TT 



and hence originates from contributions of the entire field spectrum. Usually, mechanical frequencies are much smaller 
than the resonance frequencies of the cavity: oj <^ tt/L. In this limit, the motion is slow in the time scale corresponding 
to the time-of-flight of the photon between the plates. In this quasi static limit, the susceptibility may be replaced 
by the first terms of the expansion in powers of u (we follow the notation of Ref. ||^): 

Xafsiuj) = -Kal3 + iKpUJ + fJ.al3(^^ + 0{u!^) (11) 

The meaning of the coefficients in Eq. (nW is best understood in the time domain: 



Sfait) = - ^ [Kal3&li3it) + XaflSq/sit) + /ia;3'§/3(<) + ■■•] ■ (12) 
/3=1 

The coefficients Xaf3 represent the viscous dissipative force in vacuum, in the low-frequency limit. They vanish for 
any value of D, because there are no cavity modes available at low frequencies to contribute to Imxapi'-^) as given 
by the r.-h.-s. of Eq. (0). 

In contrast to the dissipative component of the force, the dispersive component is the combined effect of all cavity 
modes even at such low mechanical frequencies. It corresponds to the coefficients Kq/j and fiap h^ (pT]) and (pT 



which may be calculated from the low frequency expansion of (10). The coefficients = — Xa/3(0) provide the linear 
correction to the static Casimir force, when the distance between the plates is changed from a to its instantaneous 
value a — 6qi{t) + 6q2{t). Hence they do not contain any new nontrivial information, and are given by 



da da? 
where 

_ r((D + i)/2) gp + i) fiL'^-^ 

^0 - (4^)P+i)/2 ^ (13) 

is the (static) Casimir energy for a cavity length a (F and C denote the gamma and Riemann zeta functions |^ ) . 
Note that in the static case, only the relative position between the plates matters, so that xii(O) = — Xi2(0). However, 
as we show below, this relation does not hold for arbitrary values of frequency. 
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Finally, are the coefficients of the inertial forces coupling the dynamics of the two plates. Replacing Eq. 
into ([lo|), we find 



D-l 



7-D-i; 



We calculate the r.-h.-s. of (Q) in Appendix B, with the help of dimensional regularization. We find 



Mil 



and 



Ml2 = 



1 TT^/^ 

67r2 2^ 



1 7r^/2 
67r2 2^ 



p-3)(i?-l)^ 



C(2-i:')+7r^/2-C(^^;L)P 



p-3)p-l) ^ 



2-f )C(2- 



-(i? + i)r 



D + l 



D-l 



hL 



D-l 



(14) 



(15) 



(16) 



Despite of the poles of the Gamma function, the first term in r.-h.-s. of ( [l5| ) and ( |l^ ) has a finite limit at positive 
even values of D larger than 2, because of the zeros of the zeta function. Those values may be obtained more easily 
from the representation given by (B6) in Appendix B. Moreover, this term has a removable singularity at D = 1, this 
time due to the pole of the zeta function, but again the limit is finite. The Casimir energy Eq, as well as the two 
inertial coefficients have negative values for any D > Q, and in Fig. 1, we plot their absolute values per unit of 'area' 
as functions of D, and taking ?i = 1, a = 1. They have minima between D = 22.8 {^12) and D = 25.2 (-^0)7 and then 
increase exponentially for larger values of D. 




In Table 1, we show the explicit expressions for _D = 1, 2 and 3. The results for _D = 1 were already derived in 
Ref. , as the linear and low- frequency limit of the exact result obtained by Fulling and Davies ^ , and also as the 
perfect-reflecting limit of the linear result for mirrors with frequency-dependent reflection coefficients. For D = 1 and 
-D = 3, we have also performed an independent calculation using an alternative regularization prescription, based on 
the introduction of a high-frequency cut-off, in order to check the results. 



TABLE I. Inertial coefficients 







D = 1 


D = 2 


D = 3 


/iii/(?ii"-7«'') 

Mi2/(7iL^~V«'') 


-(3 
-(6 


+ 7r^)/(367r) 
+ 7r2)/(727r) 


-(l + 6C(3))/(967r) 
-(-l+3C(3))/(967r) 


-ttVIOSO 
-7rV2160 



Since /in 7^ —^■i2i the dynamical Casimir forces are not functions of the relative motion only, and the effect takes 
place also in the case of rigid motion of the cavity. This fact leads to a global mass correction, which is directly 
connected to the Casimir energy as we discuss in the next section. 
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III. GLOBAL INERTIAL CORRECTION AND ITS CONNECTION WITH THE CASIMIR ENERGY. 



The inertial forces given by Eq. leads to coupled equations of motion for the two plates: 

2 

maSjait) = -^^J,a^3Sip{t) + Fa, (17) 

13=1 

where Fa represents the sum of all forces acting on plate a, excluding the inertial force itself. The coefficients flap 
go to zero as when the plates are set apart, hence uia is the free-space mass of plate a. Remarkably, the inertial 
force on a given plate also depends on the acceleration of the other plate. Even in the absence of external forces, 
in general these equations cannot be reduced to uncoupled equations for center of mass and relative motion, unless 
TOi = m2 — TO, and in this case the effective mass for relative motion is (to + /in — /ii2)/2. 

In the case of global, rigid motion of the 'cavity' system, Sqi = Sq2 = Sq, we find, by adding the equations of motion 
for the two plates, MSq" = Fcxt, where F^xt is the sum of all external forces (the static Casimir components cancel, 
as well as the external stresses which enforce the rigidity by compensating the Casimir attraction), and 

M = TOi+TO2 + 2(/iii+/fi2). (18) 

Note that fin + /ii2 is not a mass correction of a given plate taken separately, since /ii2 is associated to a coupling 
between the plates as shown by Eq. (p^). However, AM — 2(/iii -I- /ii2) is in fact a global inertial mass correction 
introduced by the Casimir effect, and from Eqs. (Esh, (ISh and (ttq) we find (now re- introducing the speed of light c) 



AM ^(D + l)^. (19) 

This relation between the inertial mass correction and the Casimir energy is in agreement with the law of inertia of 
energy for a stressed rigid body. In fact, for a plane cavity (length a) moving along the normal direction, the general 
mass-energy relation reads [ p^ 

AM ={E + aF)/c^, (20) 

where E is the energy of the system contained within the cavity, and F the corresponding force, with F > denoting 
outward force [|l^ . Here E = Eq is the unperturb( 
From the power law as given by Eq. ( p^ we find 

I _ _ 

da 



outward force [|17[. Here E = Eq is the unperturbed Casimir energy, and F = /g"'' = is the static Casimir force. 

fi''=-'^^DEo/a, 



which combined with (^0|) ver ifies Eq. (|l|). 

In view of recent measurements of the Casimir force [ p^ , we evaluate the order of magnitude of the inertial 
correction, taking D — 3 and assuming that ( |l9|) also holds for the electromagnetic case. Then the global mass 
correction per unit area is 

AM TT^h 



L2 180ca3 ■ 

For a plate separation a ~ lOOnm, we have AM/L'^ ~ —2.0 x lQ~'^'^g/cTo? , a very small correction. 



IV. CONCLUSION 



We have calculated the dynamical Casimir force for two moving plates in a D-dimensional space. We have employed 
linear response theory, which allows us to compute the dynamical force from the spectrum of fluctuations of the force 
for the static configuration. In this approach, the motion is taken to be a small perturbation, and terms of second or 
higher order in the displacements of the plates (and their time derivatives) are neglected. Moreover, we have taken 
the quasi static approximation, expanding the linear susceptibility in powers of frequency. 

In this limit, the dominant dynamical contributions are inertial forces. They have been derived from dispersion 
relations, and with the help of dimensional regularization. The inertial force on a given plate has two contributions: 
one proportional to its own acceleration (coefficient —fin), and one proportional to the acceleration of the companion 
plate (coefficient —fii2)- For D = 3, we have fin = 2/ii2, but for larger dimensions the cross-acceleration force 
dominates over the self-acceleration force, as shown in Fig. 1. The two effects jointly contribute to the global mass 
correction, which in its turn is related to the Casimir energy as predicted by the law of inertia of energy for stressed 
rigid bodies. 
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APPENDIX A: DISSIPATIVE COMPONENT OF THE FORCE 

In this appendix, we compute the spectrum of fluctuations of the force, Capiuj)^ which is directly connected by 
linear response theory to Imxapii^), as shown by Eq. (^. 

According to Eq. (^), the force operators are quadratic in the field operators. Hence, the force correlation function 
may be computed from two photon matrix elements as follows : 

cMt) = iE(oi/^"'Wk,0(e,e'l/r(o)|o), (Al) 

wher e [e , e') is a two-photon state, the label e = /, n2, njj representing a given field mode. The time dependence in 



Eq. (|Al|) is of the form 

(0|i^°)W|6,e')=expHK+c^e')i](0|/^''H0)|e,e') 
for only annihilation operators contribute in this matrix element. Thus, the spectrum of fluctuations is given by 

C„M^)='r^(0|/(")(0)|6,6')(e,6'|/f (0)10)5(^-^,-^,0- (A2) 

The matrix elements are computed from the normal mode expansion as given by Eq. (^: 

(0|/?^(0)|e,e') = ^2^^1=5„^,_„,...<S„^,_„^, (A3) 

(0|/f (0)k,6')=-vr2A^!^(_l)W'^ ^ ^..^ (A4) 



We combine Eqs. (g), (A2), (A3) and replace 

D-l 

2tt 



to derive the result for Imxii(w) as given by dT]). The result for Imxi2(cj) is very similar, except for the additional 



factor —(—!)'+' present in (A4) and not in (R3) 



APPENDIX B: DERIVATION OF THE MASS CORRECTIONS 

In this appendix, we compute the series and integral in the r.-h.-s. of (|lj), which give the results for the mass 
coefficients ^ai3- Since the integrand depends on k|| only through its modulus, we may replace the multiple integral 
by an integral over fc|| — |k|||, with 

2_-(D-l)/2 

II r(i^) " " 

The terms with li — I2 are simpler, and may be solved directly in terms of the Euler Beta function [ pO[ : 

■^^-^y^ B{{D -l)/2,-D/: 

2^+2 r((i:>- 1)/2) 



kI=1 
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and ^12(^1 = h] = — = ^2]- The series in the r.-h.-s. of (Bl) converges when D < 1 to the value C(2 — D); outside 
this range the regularized result is obtained by analytic continuation in D, with the help of the well-known analytic 
extension of the zeta function. The terms with li ^ I2 are computed in the following way: we first replace into jl^ ) 
the result 



which follows from and then the resulting integrals may also be computed in terms of the Beta function: 

(D-3)/2 1 




B 



B 



D + l D 



• "2" ) '1^2 



2/15+2 



D -1 D 



fiL 



D-l 



where ^ runs over all positive values of I2, except I2 = h- In Appendix C, we regularize the sums in Eq. (B2) once 
more with the help of the analytic extension of the zeta function. The complete result reads 



Mil 



247r 



l_ n^/^T{ D/2) g^^^ _ - 1)C(2 ~D)+ A^\D + l)C(-i?)] 



2^ 



and 



= ^iT-^ [-iD - 3){D - 2){D - 1)C(2 - D) + 2it\D + l)C(-i?)] —j^. 



247r^ 



(B3) 



(B4) 



Due to the poles of the Gamma function, these results are at first-sight ill-defined for even positive values of D. 
However, these poles are compensated by the zeros of the expression within brackets, and more useful representations 
may be derived with the help of the result 



C(i-.). 



(B5) 



We employ ( [B5[ ) for the second term within brackets in (B3), whereas for the first term we employ the recurrence 
relation for the gamma function and then use (85) with s = 2 — D. We find 



I ^-(D+l)/2 



Mil 



12 



2D 



{D-3){D-2){D-l) ^ (D - 1 
D 



C{D - I) + 2{D + l)T 



D + l 



aD + i) 



hL 



D-1 



(B6) 



and a similar representation for /ii2. However, this representation is still undefined for D = 2, because of the pole of 
the zeta function. Hence a third representation is necessary to compute this particu lar c ase. Rather than using (B5) 
for the first term, we employ the recurrence relation twice to derive, again from Eq. (B3), the representation given by 



APPENDIX C: SERIES FOR THE EVALUATION OF ^in AND /iia 



In this appendix, we compute the series appearing in Eq. (B2). Our starting point is the result 



E 

1——Q0 



1 



sin^(7ra;) tan(7rx) \ sin^(7ra;) 



2tt^x 



2^2 



(CI) 



We evaluate the Laurent series expansions of the r.-h.-s. of (CI) in powers of x — n, with n integer, and from the 
lowest order terms derive the result 
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where the star denotes that the sum is taken from 1 to infinity, the value I2 — h excluded. In a similar way we derive 

^^W^^~mt"ls- ^^^^ 

For the derivation of /ii2, we need 

l2 = l 

For the first term we change the variable of sum to fc = Z2/2 to find 

* i2 ^ * p 

^^^W^^^^^WJW^W' ^^^^ 



For even values of h this sum is given by the r.-h.-s. of (|C2| ) (where we replace h by ^i/2), whereas for odd values 
the result may be obtained directly from (jpl\ ), replacing x by li/2 (in this case, the function in its r.-h.-s. is finite at 
this value, and the sum in (|C5| ) runs over all positive integers). The final result reads 

(C6) 



Following the same approach, we compute 



After replacing these results in the r.-h.-s. of Eq. (B2), we still have sums over li to compute. The sums are divergent 
for the values of D of physical interest, but regularized expressions are obtained by taking the analytic extension of 
the zeta function. 
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